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The fast rates of magnetic reconnection found in both nature and experiments are important to understand 
theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong 
guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with 
recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, 
where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a 
more detailed derivation of the analytic model presented in a recent letter, and present additional simulation 
results to support the conclusions that the magnetic reconnection rate in this regime is independent of both 
collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and 
kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current 
layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an 
alternative interpretation for the numerical results presented therein. 
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I. INTRODUCTION 

Magnetic reconnection 1 2 is the changing of magnetic 
field-line connectivity by localised magnetic flux unfreez¬ 
ing. This process can release vast amounts of stored mag¬ 
netic energy, resulting in rapid heating and particle ac¬ 
celeration within solar flarescoronal mass ejections, 
and substorms in the Earth’s magnetosphere. 4 

Reconnection is a common process in magnetically con¬ 
fined laboratory plasmas, where it can be utilised to reach 
desirable equilibrium pressure and magnetic field pro¬ 
files 5 8 during plasma start-up. However, reconnection 
can also be undesirable in tokamaks and other devices, 
where it can destroy nested flux-surfaces and degrade 
confinement PHH 

In these applications, the details of the reconnection 
vary, depending on the magnetic field geometry and col- 
lisionality of the plasma. Low-/? reconnection, where the 
plasma-/? is the ratio of thermal to magnetic pressures, 
occurs when the in-plane reconnecting magnetic held is 
dominated by a strong out-of-plane guide magnetic held 
component. This reconnection regime is pertinent to 
tokamaks, the solar corona, and magnetically dominated 
astrophysical environments. However, despite the wide 
range and importance of these applications, reconnection 
in this regime remains poorly understood. 

An outstanding theoretical question concerns the fast 
timescales of magnetic energy release, compared with 
collisional timescales. It has been demonstrated nu¬ 
merically 12 15 and analytically®®^ that zero guide-held 
reconnection can be fast; a signihcant fraction of the 
Alfvenic rate and independent of both dissipatiorP- 4 19 20 
and system-sizep 1 22 provided that the dissipation re¬ 
gion (DR) thickness S falls below the ion-skin depth 


^Electronic mail: stanier@lanl.gov 


di = c/u p i. Here, c is the speed-of-light, uo p i is the ion 
plasma frequency, and the DR is the region where col¬ 
lisional dissipation is the dominant contribution to the 
parallel electric held associated with reconnection. With 
a hnite guide-held, it has been shown that the ion-sound 
Larmor radius p s = T e /rrii/Ul c i , where T e is the elec¬ 
tron temperature, rrii the ion mass, and Ll c i the ion cy¬ 
clotron frequency, can play the role of the ion-skin depth 
as the threshold for fast reconnection. 23 HH1 it is also worth 
noting that time-dependent studies with a guide-held 
have found faster than exponential tearing growth-rates 
in the non-linear regime when hnite electron inertiap 9 or 
hnite ion gyro-radius GHBEH effects are included. 

If the two-huid scales di or p s are large enough, the 
plasma can support fast-dispersive waves (FDWs) with 
frequency uj oc k 2 for wavenumber k , namely whistler 
and kinetic Alfven waves, respectively. It was proposed in 
Ref.|23that these FDWs play a critical role in facilitating 
fast reconnection, where the outhow velocity from the DR 
scales as the phase-velocity of the wave, v Q oc k oc 1/S, 
to give dissipation-independent electron hux through the 
DR: v Q S oc <5°. 

In Fig. 3a of Ref. |27l it was further demonstrated that 
the reconnection rate decreases with the plasma-/? for 
hxed mass-ratio and ratio of guide-to-reconnecting held. 
For very low-/?, such that /?/2 < m e /rrii with m e the 
electron mass, p s falls below the electron skin-depth, 
d e = ^ /m e /rriidi , and there are no FDWs supported 
in the plasma. Due to the decreased rate and open¬ 
ing angle between the magnetic separators, it was con¬ 
cluded that reconnection is slow in this regime. How¬ 
ever, we note that no formal scaling study against dissi¬ 
pation or system-size was presented in Ref. [27j and re¬ 
cent results have called this latter conclusion into ques¬ 
tion. Firstly, it has been demonstrated that reconnec¬ 
tion is fast in electron-positron plasmas p 3 ® which do 
not support FDWs. Secondly, kinetic simulations 37 ^ 
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have shown that reconnection remains fast in low-/3 ion- 
electron plasmas when f3/2 < m e /rrij , appearing to con¬ 
tradict earlier two-fluid results!^ 

In this paper, we address this contradiction with ana¬ 
lytical and numerical results, demonstrating that recon¬ 
nection is indeed fast within a low ~/3 two-fluid model in 
the absence of FDWs. A brief summary of our main re¬ 
sults in support of this conclusion has been presented in a 
recent letter. 39 In this paper, we give a more complete de¬ 
scription of the numerical and analytical methods, and 
present further results that support the conclusions of 
dissipation and system-size independence that we were 
unable to present in Ref. [39]due to space limitations. We 
also present new results showing detailed comparisons of 
the reconnection rate, magnetic flux, and current den¬ 
sity, between fluid and fully kinetic simulations. Good 
agreement is found between fluid and kinetic simulations 
in cases with and without FDWs, supporting the conclu¬ 
sion that the reconnection rate is independent of the DR 
physics in both cases. Finally, we perform simulations 
that can be more easily compared with those presented 
in Ref. 27] and suggest an alternative interpretation for 
the observed decrease in reconnection rate with plasma- 
/3. 


The layout of the paper is as follows. In Sec. [TTJ we 
describe the low-/? two-fluid model that forms the basis 
of the discrete analysis and numerical simulations in this 
paper. In Sec. m we discretise this model at the DR, 
and show analytically how the DR must self-adjust to 
permit fast, dissipation independent reconnection. The 
numerical method used to solve the model equations of 
Sec. [IT] is described in Sec. [IV| then in Sec. [V] we demon¬ 
strate that the DR does have the capacity to self-adjust 
in the manner required to give dissipation independent 
rates in cases with and without FDWs. In this section, we 
also compare and find excellent agreement between fluid 
and kinetic descriptions. It is further shown, in Sec. [Vl| 
that rates are a significant fraction of the Alfvenic rate 
and independent of system-size, so that reconnection is 
formally fast both with and without FDWs. The compar¬ 
ison with Ref. EH is given in Sec. |VII[ before conclusions 
are drawn in Sec. Enin 


II. LOW-/3 TWO-FLUID MODEL 

The low-/3 (yd<C 1), two-dimensional (c^=0, where 
z='Vz is the o ut-of-pla ne direction), two-field reconnec¬ 
tion modeP ^ 26 E 8 l 4Q l 41 J can be derived from the two-fluid 
magnetised plasma equations^ assuming uniform tem¬ 
peratures with cold-ions (VT^=VT e =0 and T^<CT e for 
ion temperature T^), strong-guide field (|£?|<cF>o where 
B=zx'V r i/j is the in-plane field, -0 is the flux, and Bqz is 
the guide-field), an MHD ordering (the electric drift ve¬ 
locity is on the order of the ion thermal speed Ve VTi ) 
and small ion parallel flow (v\\i<^vri)- As is typically 
done (see e.g. Ref. [43]) . we also assume a simple per¬ 
pendicular viscosity closure for both ions and electrons. 


Here, the equations are given in terms of the vector mag¬ 
netic field B for the analysis that follows, and are nor¬ 
malised by a characteristic Alfven velocity and macro¬ 
scopic length-scale. 

(d t + v ■ V)w = B ■ Vj + /zV 2 w, (1) 

d t B* - V x (v x B*) = p 2 s V x [B x (z x Vw)] (2) 

— V x [V x ( r/B - r/H^7 2 B)] , 

where v = z x V</> is the in plane velocity in terms of 
streamfunction </>, u = V 2 </> is the vorticity, j = V 2 ^ 
is the parallel current density and B* = £? + d 2 V x 
(V x B). The normalised dissipation coefficients are the 
ion collisional viscosity /i, the plasma resistivity r/, and 
the hyper-resistivity t]h that can be written in terms of 
an electron collisional viscosity p e as t]h = d 2 p e . 

When d e =r]=r]H=0, Eq. © can be written in flux 
form as dt'ip+v s -'V'ip=0 1 where v s =zx'V ((/)— p 2 V 2 0) is 
the velocity that advects the frozen-in magnetic flux. 26 44 
Physically, this is an electron perpendicular velocity that 
combines the lowest-order electric drift v (the single-fluid 
velocity in this model) and an electron diamagnetic drift 
that can become significant close to the X-p oint due to 
the characteristic density asymmetry patternPEHl associ¬ 
ated with strong guide-field reconnection. 

Relaxing the d e =0 condition gives dt'ip-\-v s -'V r ip = 
d 2 e ( d t j + v • Vj). Here, in the last term on the right, 
we have v • Vj rather than v s • Vj, since the diamag¬ 
netic contribution to electron inertia cancels out with the 
lowest order electron gyro-viscosity (see Ref. [45]) . Thus, 
we avoid the potentially large and spurious term that 
is present in non-reduced Hall-MHD simulations that 
include electron inertia while neglecting electron gyro- 
viscosity. 

In a uniform (B = B 0 z) collisionless (r]=fi=r]H=0) 
plasma, Eqs. 0 and 0 support waves with dispersion 
relation 


W = fc|| V / (l + p2fc 2 )/( 1 + < i2 fc 2) ; (3) 

where k\\=z-k. In the limit p 2 s k 2 1 d 2 & 2 , Eq. ^ 

describes the fast-dispersive kinetic Alfven wave with 
uj=p s k\\k. However, for p s < d e there are no FDWs. 
There is only one condition here for the presence of 
FDWs, in place of the two-conditions discussed in Ref. 27 ; 
since the condition for the presence of whistler waves 
has been ordered out due to the strong guide field as¬ 
sumption. The condition p s < d e is equivalent to 
/d/2 < m e /mi. 

III. DISCRETE MODEL 

A summary of this section can be found in the recently 
published letter. 39 ! Here, we additionally include the in¬ 
termediate steps in the derivation of the governing equa¬ 
tion for the DR thickness, and highlight the similarities 
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FIG. 1. Dissipation region (DR) embedded in a larger re¬ 
gion of thickness h = max[p s ,d e ] at which two-fluid effects 
become important. The length w is defined as the extent 
of the v sx outflow jets, and w v the extent of the v x outflow. 
Reprinted with permission from Stanier et al., Phys. Plasmas, 
22, 010701, (2015). Copyright (2015) American Institute of 
Physics. 


Two-fluid region 



between strong guide field reconnection in ion-electron 
plasmas without FDWs, and reconnection in the low ~/3 
pair-plasma regime studied previously in Ref. [35] 

To construct the analytic model, Eqs. 0 and 0 are 
discretised at the DR as shown in Fig. [l] using the tech¬ 
nique of Refs. Mi EHI andl35l The discrete magnetic field 
components are defined on the edges of a rectangular DR 
of thickness S and length w, so that B x = x • £?(0,£/2), 
By = y • B(w/ 2,0), and the discrete flow stream func¬ 
tion is <F = —<fi(w/2,5/2). Then, the inflow and outflow 
velocities are given by V x = 2<F/£, V y = —2 $/w. It 
is assumed that the DR is quasi-steady, so that deriva¬ 
tives with respect to time are small compared with other 
terms and can be neglected. This is normally true when 
the system reaches non-linear saturation, such as at the 
time of peak reconnection rate in the island-coalescence 
problem as shown in Sec. [V] and Ref.[46| Also, neglecting 
numerical factors of order unity gives three equations for 
the five unknowns (S, w, B x , B yi <F) 



Defining £ = S/w, b = B y /B x , p s = p s /S and d e = 
d e /S, Eqs. ( [5|6| ) give 

b = 1 + P 2 S (l + £ 2 ) + 2d 2 

i + p2(i + £ 2 ) + 2<& 2 ’ 

which, after substitution back into Eq. ([6|, gives an equa¬ 
tion for the discrete stream function 

* = D ((-' -() 

i + (i,i + 4)(i + ( 2 ) 

Without loss of generality, we consider <f> > 0 so that 
£ < 1. Substituting Eqs. 0 and 0 into Eq. 0 gives 
the master equation for £ 


[l+tf (I+O] 

i +1 

[p 2 s+dl) 

(i + of 


\ + pi(i + ?) + 2die 

1 2 


(i-G 2 

S 2 £ 4 


(l + a 2 [l + (# + <£) (l + O 


SS^ 


(9) 


where S = \[2B x w/D is the Lundquist number and S y = 
y/2B x w/p. This equation can be solved numerically for 
£ (and therefore (5) for given parameters B x and w. The 
reconnection rate (E z ) in steady-state (d t j = 0) can then 
be calculated using the discrete expression 


E a D (*k - By) * (r 1 - e) 

V w ) ~ l + 2d 2 £ 2 / [1 + p 2 a (1 + e 2 )] ’ 

( 10 ) 

Here we solve for £ and E z analytically, noting that large 
E z preferentially occur for £ 1. This motivates the 

approximation l + £ 2 ~ 1 — £ 2 ~ 1, which considerably 
simplifies the above expressions. We also assume that 
2d 2 £ 2 = 2 d 2 e /w 2 < 1. The validity of these two assump¬ 
tions was checked for the numerical simulations presented 
in this paper. Finally, in this paper, we concentrate on 
viscous DRs (rj = 0; t]h, H 7^ 0), as resistivity alone can not 
prevent the DR collaps ing to zero thickness once S < h, 
where h = max[p s , d e ]!^^ 7 The resulting equation for 
S{B X , w) is given by 


£ 4 ((5 2 + ^ + Q 3 

w 8 S 2 + p 2 


SI 


+p 2 +d 2 ) 

rjH 


( 11 ) 


where Sh = V%B x w 3 /tjh is the hyper-resistive Lundquist 
number, and the quasi-steady (dtj = 0) reconnection rate 
is 

E x * D (A - « V hB x /S 3 . (12) 


where D = p + t]h A, and A = S~ 2 + w~ 2 . This set of dis¬ 
crete equations combine the finite-p s terms in Ref. [28l 
and the finite-d e terms from the discrete equations in 
Ref. m Eqs. @-(§) are invariant under plasma flow 
reversal about the DR (B x , B y ,Q,5) (B y , B x , — <L, w). 


A. Single fluid scalings 

We fi rst note that the usual resistive single-fluid 
scalings SHED including the correction for finite ion- 
viscosityj^ can be straightforwardly obtained from 
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Eqs. ([9 |[To| as was shown in Ref. |28j Concentrating here 
on hyper-resistive layers with Eqs. cco and taking 
the single-fluid limit (S 2 d 2 ,p 2 ) with zero ion-viscosity 

(m = 0 ), gives 

S = 5 H = w Sh X/ \ E z = E zH = V2B 2 x Sh 1; \ (13) 

which are the expected single-fluid hyper-resistive scal¬ 
ings for the DR thickness and reconnection rate (see e.g. 
Ref. f5T|) . Here, as p s = 0 , there is no separation between 
in-plane ion and electron flows, v = v s , so that the DR 
length w is equal to the extent of the single-fluid outflow 
jets, w v , see Fig. [lj These scalings can also be generalised 
for finite ion-viscosity . 39 


B. Scalings for case without FDWs 


In what follows, we assume for simplicity that the ion- 
viscosity is set to a fixed value of p/ph = l/(Ps + d 2 ). 
However, the following two-fluid scalings can be straight¬ 
forwardly generalised to relax this assumption. 

We consider first the two-fluid case without FDWs, 
with d 2 p 2 ,S 2 . In the limit of small-p s , such that 
d 2 ^5 2 ^> p 2 , Eqs. (11 12) give the DR thickness and 
reconnection rate as 


(5 = 


Phw 

B x dV 


E z = 


pjjW 3 


(14) 


The reconnection rate appears to be “super-fast”, that 
is E z oc t]h with a < 0. However, the free parameters B x 
and w have not yet been specified, and they do scale with 
Ph as will be described. 


The expression for the DR thickness in Eq. (14) is 
similar to that for a low ~/3 electron-positron plasma, 
S=pw/(B x d e ), in Refl35 Indeed, for the electron- 
positron mass ratio, rrii=m e , such that pn=d 2 p e =d 2 p , 
we obtain the elec tron -positron result. This is because 
the two-field Eqs. © reduce, in the limit rrii=m e and 
p s = 0 , to the low -f3 electron-positron plasma equations 
studied in Ref. |35j apart from numerical factors of order 
unity. 

As discussed in Ref. EH 2 the electron viscous (hyper- 
resistive) DR supports the reconnection electric held at 
the X-point in steady-state, but the bulk of the cur¬ 
rent is supported at d e -scale such that the magnetic 
held at the DR edge scales like B x « 5B x d/d e , where 
B x d = x B( 0, d e f 2) is shown in Fig. [l] Finally, as p s is 
small in this limit, there is little separation of ion and 
electron in-plane hows, and so w ~ w v . This gives the 
following scalings for the DR thickness and reconnection 
rate 



E z 


d e 

'Wv 


TD 2 

n xd~r ? 


(15) 


where 5d = w v (ph / B x dW 3 ) x / 4 is equal to the single-huid 
thickness in Eq. (fl3[) , evaluated with B x d and w v . E z has 


no explicit dependence on hyper-resistivity in Eq. (15), 


but it is still a function of the parameters B x d and w v 
which may have some pn dependence. It is shown in 
Section [V] that this rate is indeed dissipation independent 
as these parameters do not depend on 77 #, in agreement 
with the electron-positron result.^ 

In the second limit d 2 >> p 2 S 2 , Eqs. (11, 12) give 
the DR thickness and reconnection rate as 


( wtjhPs \ 1/2 

Miwr) 


E x = 


B*4 

PhPs w3 


1/2 


(16) 


The rate again would be super-fast if the free-parameters 
B x ,w oc p® H . However, the rate can not become singu¬ 
lar as ph 0 , so we expect it to remain dissipation- 
independent when finite-p s effects become important, 
provided that the DR has the capacity to self-adjust B x 
and w appropriately. In this limit, the current density is 
still supported at d e -scale, so B x scales as before. Com¬ 
paring Eqs. ([16, 15), the DR length w must adjust as 
w = ( S/p s )w v to keep the same rate. In Section [V| we 
verify that the DR does have the capacity to self-adjust in 
this manner, even in the absence of fast-dispersive waves, 
to give the DR thickness and dissipation-independent 


rate of Eq. (15) that is constant across both limits. 


C. Scalings for case with FDWs 


Next we consider the two-fluid case with FDWs in two 
limits; p 2 d 2 ^> 5 2 , and p 2 S 2 d 2 . In both limits, 
the DR thickness and reconnection rate are found from 


Eqs. (11, 12) as 


fr/ H w\ 1/2 


E z = 


(B*p\ 


5J\l/2 


\Ph 


(17) 


where again E z appears proportional to a negative power 
of ph • Previous studies^^ have shown that reconnec¬ 
tion is indeed dissipation independent in this regime, and 
so we expect that the DR is able to self-adjust to give a 
fast rate, as in the d e > p s case. In the first limit 5 < d e 
so that B x oc 5B x d/d e as before, and a rate that is ex¬ 
plicitly independent of pn requires that w oc S. We show 
from numerical simulations, in Sec. [Vj that the fixed- 
aspect ratio scaling S/w oc p s /w v , from the d 2 ^> p 2 ^> S 2 
limit, holds whenever S < p s , regardless of the value of 

d e • 

The second limit has 5 > d e , and so a scaling for B x is 
not known a priori from physical arguments. However, 
with the w = ( 5/p s )w v scaling for S < a dissipation 
independent rate in this limit also requires that B x oc S. 
From simulations in Sec. [V] we find for both limits that, 
to a good approximation, B x ~ ( S/p s )B xs with B xs = 
x • £?(0,p s /2). With this, the DR thickness and rate in 
both limits is given by 


S 



T7i _ jD 2 rS 

Ez — £> xs , 


w v 


(18) 
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where S s = w v (t]h /H^u; 3 ) 1 / 4 . Again, this rate has no 
explicit dependence on 77 #, and we verify in Sec. |V|that 
the two-fluid region parameters B xs and w v have no t]h 
dependence. 

The DR thickness and reconnection rate in cases with 
and without FDWs can be expressed simply as 

^ = ~j^ 5 E z = v y hB x h = B xh ——, (19) 


where v y h = B x hh/w v is the inflow velocity into the 
larger two-fluid region, see Fig. [l] constrained by flow 
continuity with the outflow equal to the upstream Alfven 
speed ( B x h in these normalised units). This form of the 
reconnection rate can be arrived at from simple steady- 
state arguments, considering only the outer two-fluid re¬ 
gion. But, by starting from a discrete model of the DR, 
we have illustrated how S, B x and w must scale in each 
case to give this result. 


The form of the rate in Eq. (19) also holds for the 


single-fluid limit, when h = 5h and w = w v , sim¬ 
ply because the reconnection rate must be the same 
as the upstream rate of inflowing flux in steady-state 
(d y E z = —dfB x = 0). However, the mechanism whereby 
this result is obtained differs between single and two-fluid 
regimes. In the single-fluid regime, the reconnection rate 
is limited by the DR, so that the upstream rate of in¬ 
flowing flux adjusts and slows down to the DR rate in 
steady-state. This is in contrast to the two-fluid regime, 
where the rate is instead bounded by rate of flux inflow 
to the outer two-fluid region, and the DR self-adjusts to 
match this upstream rate. The crucial difference gives a 
rate explicitly dependent on t]h in the single-fluid case, 
and a rate independent of DR physics in the two-fluid 


strict control of numerical dissipation in the scheme used, 
and no additional high-order dissipation than is stated 
explicitly in Eqs. (|TJ |2j, we are able to accurately mea¬ 
sure how DR parameters such as the thickness, length 
and upstream magnetic field depend on the magnitude 
of the physical dissipation. Both grid and timestep con¬ 
vergence studies were performed for a selection of sim¬ 
ulations to verify the accuracy of the results presented. 
Finally, we note that the solver successfully reproduces 
the linear tearing eigenmode structure and growth-ratd^ 
in the two-fluid (large p s regime), as calculated by a sep¬ 
arate linear stability code. 

The simulations presented in this paper use one of two 
typical reconnection problem set-ups; the island coales¬ 
cence problem, and the single Harris-sheet. The island 
coalescence simulations have initial flux-function 

'ipic = A In [cosh (x/X) + e cos (y/ A)], (20) 

with A = l/(27r), e = 0.2, and with a sinusoidal pertur¬ 
bation of 0.1%. This problem is solved in a quarter box 
(x,y) G [0,1] x [0,1] with symmetry (or anti-symmetry, 
depending on the variable) boundaries at x = 0, y = 0 
and y — 1, and a perfectly conducting boundary at x = 1. 

The Harris-sheet runs use the flux-function 

Ha = -A In [cosh (y/ A)]. (21) 

They are run in a half-box (x,y) G [—L X ,L X ] x [0,0.5], 
with equilibrium current sheet thickness A = 1/(87r), and 
initiated with a sinusoidal perturbation with magnitude 
equal to 3% of the upstream boundary magnetic field 
strength. The half-box domain ensures that any sec¬ 
ondary islands that form at x = 0 are not trapped in¬ 
definitely. 


IV. NUMERICAL METHOD 

We solve implicitly the vorticity and magnetic-flux for¬ 
mulation of Eqs. with a Ja^obian-Free Newton 

Krylov (JFNK) based algorithm To accelerate con¬ 
vergence, a physics-based preconditioning strategy was 
developed, where the hyperbolic couplings present due 
to the fast-dispersive wave terms in the preconditioner 
are parab olise d before they are inverted using multigrid 
methods! 53 * 54 * This preconditioning has no effect on the 
solution, but reduces the numerical stiffness when taking 
larger time steps than the explicit step associated with 
the fast-dispersive waves. The result is a large effective 
speed-up over explicit schemes. 

The spatial operators in Eqs. are discretised us¬ 

ing second-order centred finite differences, apart from the 
advection terms, which are discretised using quadratic 
upstream interpolation (QUICK). For timestepping, we 
use either second-order Crank-Nicolson or Backward- 
Differentiation formula (BDF-2). We use grid-packing 
focused around the X-point to adequately resolve the 
DR, which is often sub-d e scale, in all simulations. With 


V. DISSIPATION INDEPENDENCE 
A. Reconnection rates 

Figure shows the flux and current profiles in the 
initial conditions (t = 0 , left), and at the time of peak 
reconnection rate (t = 6 . 8 , right) for an island coales¬ 
cence run with d e = 0.01, p s = 0.002 (without FDWs), 
tjh = 10 -9 and p = 10 -5 . At the time of peak rate, 
the X-point has collapsed and a current-sheet (blue) has 
formed between the two islands. The reconnection rate 
is plotted against time in Fig. [2 ]d, where it has been 
normalised by the Alfvenic rate at the upstream bound¬ 
ary (x = 1). The initially ideal island coalescence in¬ 
stability quickly reaches non-linear saturation when the 
rate peaks, at which time-dependent terms such as the 
electron inertial contribution to the X-point electric field 
(blue) are typically small. For this reason, we use island 
coalescence simulations to benchmark the quasi-steady 
discrete model of the DR presented in Section m 

Figure [ 3 ] shows the peak rate from a number of sim¬ 
ulations that are identical to that shown in Fig. [2j but 
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FIG. 2. a) Current density (colour scale) and magnetic 
flux contours (black dotted) for an island coalescence simu¬ 
lation with d e — 0.01, p s — 0.002, t\h — 10 -9 , p — 10 -5 at 
t — 0 (left) and at peak rate (t — 6.8, right). The quarter- 
domain simulation has been reflected in x — 0 and y — 0, 
and then cropped at \x\ = 0.5. b) Reconnection rate (E z , 
black) against time for this simulation, with inertial (blue) 
and hyper-resistive (red) contributions. The vertical dashed 
line indicates the time of peak rate. 



FIG. 3. Peak reconnection rate against hyper-resistivity from 
island coalescence simulations with d e = 0.01, p s — 0.002 and 
p/p h = 10 4 . The relative size of the DR thickness 5 with 
respect to d e and p s at the time of peak rate is indicated. 



FIG. 4. Peak rates from island coalescence runs with d e = 
5 p s = 0.01 (orange A), d e = 0.01, p s = 0 (green x), p s — 
5 d e — 0.01 (blue O)? Ps — 0.01, d e — 0 (red +), and p s — 
d e — 0 (black y). Reprinted with permission from Stanier 
et al., Phys. Plasmas, 22, 010701, (2015). Copyright (2015) 
American Institute of Physics. 


with different hyper-resistivity and ion-viscosity. For this 
scaling study, we keep the ratio of ion-viscosity to hyper¬ 
resistivity fixed as p/rjH = l/(Pg + d 2 ) ~ 10 4 to com¬ 
pare with the scalings presented in Sec. m As hyper¬ 
resistivity is reduced between simulations, the DR thick¬ 
ness 6, measured at the peak rate, decreases so that it 
falls below d e -scale for tjh < 10 -8 and also falls below p s 
for tjh < 10 -9,5 . It is clear for this case without FDWs 
that the reconnection rate is independent of the dissipa¬ 
tion, t}h , provided that 5 < d e . In addition, there is no 
change when S falls below the smaller of the two-fluid 


scales, p s , as anticipated in Sec. IIIB 


Figure[4]also shows the peak reconnection rates against 
hyper-resistivity from a number of island coalescence sim¬ 
ulations. Shown are two-fluid cases without FDWs (or¬ 
ange A for d e = 0.01 and p s = 0.002, which are the same 
as in Fig.[3j and green x for d e = 0.01 and p s =0), with 
FDWs (blue O f° r Ps = 0-01 and d e = 0.002, and red + 
for ps = 0.01 and d e = 0), and single-fluid runs (black 
y, where p s = d e = 0). The single-fluid rate depends 
on dissipation as E z oc p^ 25 , which is given in Eq. (13) 
with B x ,w oc rf H in the single-fluid regime. However, all 
two-fluid runs with and without FDWs are independent 
of tjh- Also, there is no appreciable change in the rate 
between d e = 0.01 simulations (orange A and green x) 
and ps = 0.01 simulations (blue O an d red +), and the 
rate is also unchanged when the smaller of the two-fluid 
scales is set to zero (green x and red +). This is in agree¬ 
ment with the scaling of Eq. (19), where only the larger 
of the two-fluid scales h = max[ps,d e ] is explicit in the 
reconnection rate. 


All simulations discussed so far have used the island 
coalescence problem set-up. Fig. [5] shows the reconnec¬ 
tion rate against time for Harris sheet simulations with 
d e = 0.01, p s = 0.002 for different p#. Here, the rate is 
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FIG. 5. Normalised reconnection rate E* vs. time for Harris- 
sheet simulations with d e = 0.01, p s = 0.002 for different 
hyper-resistivity. 


measured as 

F* = 1 Wr (99] 

z (x ■ v A )(x ■ B) dt ' y ’ 

where = max( / 0 ) — min^) is the flux difference be¬ 
tween the X and O-points, and (x • va)(x • B ) is the 
Alfvenic rate at 4 d e upstream of the main X-point. In 
the series of simulations shown in Fig. [5j there is no sec¬ 
ondary island formation. The reconnection rate is again 
independent of 77 #, suggesting that dissipation indepen¬ 
dence is a general property of reconnection in the two- 
fluid regime without FDWs, rather than being specific to 
the island coalescence problem set-up. 



0.0 0.2 0.4 0.6 0.8 1.0 


6 /h 

FIG. 6. Scaling of B x (hollow markers with linear fits) and 
B x h (small filled markers) against S/h, where h — max[p s , d e ], 
for d e ss= 5 p s = 0.01 (orange A) and p s = 5 d e = 0.01 (blue O)- 
Reprinted with permission from Stanier et al., Phys. Plasmas, 
22, 010701, (2015). Copyright (2015) American Institute of 
Physics. 



B. Discrete model scalings 


The scalings for B x and w given in Sec III 


are now 


verified using the previously described island coales¬ 
cence simulations. For each simulation, we measure 
(5, w,B x ,w v , B x h) at peak reconnection rate. The thick¬ 
ness, 5, is measured in the same way as in Refs. I35landl39l 
which is appropriate for viscous DRs. The lengths w and 
w v , shown in Fig. [TJ are defined as the distance between 
the maxima of \x • v s \ y= o and \x • v\ y= o , respectively. 
The magnetic field strengths B x and B x h are evaluated 
at (x,y) = (0,5/2) and (x,y) = ( 0 , ft/ 2 ) respectively, as 
shown in Fig. []] 

In the discrete model scalings of Sec. m the scaling 
B x oc ( 5/d e )B x d was used for S < d e , as the current is 
supported at d e -scale and magnetic field can not pile-up 
on scales below this. Also, when p 2 s S 2 d^, it was 
required that B x oc S to give a dissipation-independent 
rate. Fig. [ 6 ] shows how B x , and the magnetic field at the 
edge of the larger two-fluid region B x h (see Fig. |T]) scale 
with S/h, where h = ma x[p s ,d e \. The B x parameter 
decreases with 5/h, compared with B x h that is approxi¬ 
mately constant over the full range of S/h. Furthermore, 
the trends of B x and B x h and the absolute values are 
similar in cases with and without FDWs. Here, the DR 


FIG. 7. Scaling of w (larger hollow marker) and w v (smaller 
filled marker) against hyper-resistivity for runs with d e = 
5 p s = 0.01 (orange A). The relative size of the DR thick¬ 
ness S with respect to d e and p s at the time of peak rate is 
indicated. 


parameter B x has been fit with two straight lines for both 
cases, with and without FDWs, as there is a change in 
slope by ~ 2 when S falls below the smaller of the two- 
fluid scales. This change in slope is not captured by the 
discrete model described in Sec. El as numerical factors 
of order unity are not retained in the discretisation of the 
two-field equations. 

In the limit d 2 p 2 S 2 of the case without FDWs, 
it was argued that the DR length w must scale with con¬ 
stant aspect-ratio as re oc ( S/p s )w v when S < p s , to give 
a constant and dissipation independent rate across the 
limits d 2 > S 2 pi and d? e p 2 s S 2 . Fig. 0 shows 
w and w v at the time of peak reconnection rate plotted 
against t]h for the island coalescence simulations with 
d e = 0.01 and p s = 0.002 (the case without FDWs that 
is shown in Fig. [3|. For p s < S < d e the DR length 
w is approximately w « w v , as might be expected from 
the definition of the flux-carrying velocity, v s . However, 
for S < p s < d e , where p s effects become important, w 
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FIG. 8. Scaling of w against (5/p s )w v from d e = 5 p s = 0.01 
runs with 5 < p s (orange A), and from all runs with p s = 
5 d e — 0.01 (blue Q)- 



FIG. 9. Ratio of measured DR thickness 5 to that from dis¬ 
crete model in Eq. ( fl9| for island coalescence simulations with 
d e = 5 p s = 0.01 (orange A) and p s = 5 d e = 0.01 (blue O)* 
Reprinted with permission from Stanier et al., Phys. Plasmas, 
22, 010701, (2015). Copyright (2015) American Institute of 
Physics. 


begins to shrink with respect to w v . Across both limits 
w v is independent of tjh , as is required for a dissipation 
independent reconnection rate in Eq. (15). 

Figure [ 8 ] shows the values of w from the simulations 
with 5 < p s < d e in Fig. [t] plotted against (S/p s )w v . 
There is a clear linear scaling of w with S (quadratic in tjh 
for Fig. [7]), which verifies that the DR self-adjusts in the 
manner required to maintain a constant and dissipation 
independent rate betw een th e limits d? e S 2 >> p\ and 


_ . 

4 > Pi » <5 2 (see Sec. |lIIB[ ). Also shown in Fig. [ 8 ] is the 
constant aspect-ratio DR scaling of w oc ( S/p s )w v from 
the case wit h FD Ws (p s = 0.01, d e = 0 . 002 ) that was 
Here, there is a small change (~ y/2) 


me 


used in Sec. _ 

in the slope of the linear fit between the 4 » P 2 S (no 
FDWs) and pi d 2 e (FDWs) cases that is not captured 
in our simple analytic model, as we lose numerical factors 
of order unity in the discretisation process. 

Finally, Fig. [9] shows a comparison of the measured 
DR thickness S with the theoretical scaling of Eq. (19) 
for the cases with and without FDWs. The theoretical 
and measured values agree within a factor of ~ 1.5 over 
almost four orders-of-magnitude in 77 #, or almost two 



FIG. 10. Normalised reconnection rate E* against time for 
Harris-sheet run without FDWs ( d e = 5 p s = 0.01), from a 
fluid simulation (black) and kinetic simulation (red). The 
fluid simulation is the 77 h — 10 -9 run from Fig. [ 5 J 


orders-of-magnitude in 5 (as S oc 77 J/ 2 ). 

As was discussed in Sec. mm the self-adjustment of 
the DR appears necessary to prevent unphysical “super¬ 
fast” rates. Here, we have verified that the DR self¬ 
adjusts in the manner required to match the rate of flux 
inflow to the larger two-fluid region, which is an upper 
bound on the possible rate, rather than adjusting to give 
a slow (DR limited) rate. 


C. Comparison with kinetics 

It has been demonstrated in cases with and without 
FDWs that the DR has the capacity to self-adjust, such 
that the reconnection rate becomes independent of DR 
physics. Therefore, reasonable agreement might be ex¬ 
pected between fluid and kinetic rates within the regime 
of validity for Eqs. ([lj 2L namely with cold ions and 
strong guide-field. Fig. 10] shows the reconnection rate 
from the L x = 0.5 Harris-sheet simulation of Fig. [5j 
with d e = 5 p s = 0.01 (no FDWs), \± = 10 -5 and 
tjh = 10 -9 , and a fully kinetic VPIC 5 ® simulation us¬ 
ing force-free initial conditions and with the same d e 
and p s . Additional VPIC specific parameters used are 
a mass-ratio rrii/m e = 15.2, guide field to reconnecting 
field ratio B g = 17.4, electron thermal speed to speed-of- 
light ratio v t h, e / c = 0.0625, the ratio of electron plasma 
frequency to gyrofrequency using the upstream field is 
(jj pe /£} C ex = 10 , and the ratio of ion to electron tempera¬ 
ture is Ti/T e = 1/4. There is remarkable agreement be¬ 
tween the VPIC and fluid rate, which demonstrates that 
the essential physics is captured in Eqs. and that 

the rate does not depend on the DR physics, which dif¬ 
fers between the two descriptions. This latter point can 
be seen in Fig. mi which shows a comparison of the cur¬ 
rent density and flux evolution between these fluid and 
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FIG. 11. Snapshots of current density (colour) and magnetic flux (contours) from fluid (PIXIE2D) and kinetic (VPIC) simu¬ 
lations at t = 0 (top), t = 1.6, t = 2.64 and t — 3.56 (bottom). The set-up is the same as for the Harris-sheet reconnection 
simulations in Fig. [5] with d e — 5 p s = 0.01. The current density is normalised by the value at the null-line at t = 0. 


kinetic runs. The kinetic runs show different electron 
scale physics, such as a kinking of the electron outflow 
jet and formation of a small secondary island at t = 1.6, 
which are not present in the fluid runs. However, at 
larger scales the two runs look very similar and, in par¬ 
ticular, the evolution of the layer length is the same in 
both descriptions. 

A number of studies have compared fluid and kinetic 
reconnection rates in the p s > d e regime. Ref. f56l showed 
reasonably good agreement (within a factor of \/2) in 
the non-linear rates between a cold-ion gyro-kinetic sim¬ 
ulation and a (non-reduced) Hall-MHD formulation with 
electron inertia. Ref.[57]also found good agreement in the 
reconnected flux over time between a fully kinetic VPIC 
simulation and a Hall-MHD formulation with electron 
inertia. Ref. [58] presented a comparison between gyroki- 
netic and reduced gyro-fluid simulations of the tearing 
instability. In the non-linear regime, excellent agreement 



Time 

FIG. 12. Normalised reconnection rate E* against time for 
Harris-sheet run with FDWs ( p s = 5 d e = 0.01), from a fluid 
simulation (black) and kinetic simulation (red). 
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FIG. 13. Snapshots of current density (colour) and magnetic flux (contours) from fluid (PIXIE2D) and kinetic (VPIC) Harris- 
sheet simulations with p s m 5 d e — 0.01 at t — 0.71 (top), t = 1.4, t = 2.11 and t = 2.8 (bottom). The current density is 
normalised by the value at the null-line at t — 0. 


was found between the two descriptions in the growth 
and saturation of the magnetic islands. Finally, Ref. [59| 
found excellent agreement between reduced MHD and a 
novel fluid-kinetic model for the magnetic island flux at 
saturation in the large A' regime. 

Figures M and [13] compare the reduced formulation 
in Eqs. 011 and fully kinetic VPIC simulations in 
the force-free Harris-sheet set-up of Fig. El but with 
p s = 5 d e = 0.01. Additional fluid parameters are t]h = 
10 -9 and p = 10 -4 , and kinetic (VPIC) parameters are 
ttt-2 /m e = 380.3, Bg = \‘2, Vfj l ^ e fc = 0.308, ujp e /£l cex = 10, 
and Ti/T e = 1/10. The plasma /3 calculated with the in¬ 
plane magnetic field is large in this kinetic simulation, 
and so 16, 000 particles per cell were required to reduce 
the noise to an acceptable level. The normalised recon¬ 
nection rate in both fluid and kinetic runs is E *, the 
same as in Eq. (22), normalised by the Alfvenic rate at 


a distance of 4 p s upstream of the dominant X-point. 


Fig. 12 shows at early time (t ~ 0.6) that the kinetic 
rate has a slightly sharper increase than the fluid rate, 
and at late time the reconnection saturates slightly earlier 
in the kinetic simulation (t = 3.3) than the fluid simula¬ 
tion (t = 3.65). However, the overall agreement is very 
good, and in particular the peak rates are in agreement 
to within 10%. 

The slight differences earlier on may be due to a sec¬ 
ondary island that forms in the kinetic run at t ~ 0.5 and 
is ejected at t ~ 0.7. This island can be seen in Fig. 13 at 
t = 0.71. After this island ejection, the fluid and kinetic 
plots show good agreement in the structure of the current 
layer. Also, in both cases with and without FDWs, there 
is is reasonable agreement in the magnitude of the peak 
current density between fluid and kinetic runs. However, 
this quantity does depend on the dissipation, as we see 
larger peak current density in fluid runs with the lower 
hyper-resistivity pn = 10“ 11 (not shown). 
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It is also interesting to compare the structure of the 
current layer in Fig. 13 where d e = 5 p s = 0.01, with 
Fig. fl3j where p s = 5 d e = 0.01. In the former (without 
FDWs) the current layer is both thicker and longer than 
in the case with FDWs. Finally, we note that both the 
peak rates and the time at which reconnection saturates 
is within a factor of 2 between Fig. [10| and Fig. [l2j 


VI. SYSTEM-SIZE INDEPENDENCE 


Historically, fast-reconnection has been defined as hav¬ 
ing a rate independent of the Lundquist number,® which 
corresponds to a rate independent of both collisional dis¬ 
sipation and system-size in the single-fluid regime. How¬ 
ever, in the two-fluid regime, additional physical length 
scales are present and it is necessary to check the dissipa¬ 
tion independence and the system-size independence sep¬ 
arately. Here we present a system-size dependence study 
in the low ~/3 regime, considering the cases with FDWs 
(p 2 d 2 ) and without FDWs ( d 2 ^ p 2 s ). For this study, 
we use the Harris-sheet reconnection problem set-up, as 
it is not yet clear whether there is good agreement be¬ 
tween fluid and kinetic simulations of island coalescence, 
particularly for large islands, see e.g. Ref. [Ml 

The requirements for a system-size independent recon¬ 
nection rate can be understood by considering how the 
local two-fluid region parameters in Eq. ( [l9| ) vary with 
system-size, assuming that both the DR and outer two- 
fluid regions are approximately in steady-state at the 
time the rate is measured. For the study presented in 
this section, we vary the length of the equilibrium cur¬ 
rent layer, L x , and keep the box size in the inflow di¬ 
rection, L y , and the ratio of the two-fluid scale to equi¬ 
librium current sheet thickness, h/ A, fixed. With this, 
the upstream field B x h does not vary appreciably with 
L x , and so the condition for system-size independence 
becomes w v oc L?. If there is no mechanism to limit w v . 


then Eq. (19) suggests that the rate will decrease with 
system-size as E z oc w~ x oc L~ l . 


Figure 14 shows the normalised reconnection rate E*, 
see Eq. (22), against time from Harris-sheet fluid simu¬ 


lations with the same set-up as in Figs. [5j [lOj but with 
varying L x . In contrast to the L x = 0.5 run, where there 
is no secondary island formation, the sharp increase in 
the rate for the L x = 2 and L x = 8 runs occurs along 
with the formation of secondary islands, at t ~ 3.2 and 
t ~ 5 respectively. At later times there is significant 
secondary island formation in both of these runs (see be¬ 
low), and the rates are somewhat spiky. Although the 
runs with larger L x take longer to reach the initial sharp 
increase in the rate, there is no clear dependence of the 
rate on the system-size L x after this onset. 


Figure 15 shows the peak reconnection rates < E* > 


plotted against system-size, with E* as in Eq. (22) and 


normalised by the Alfvenic rate at 4 h upstream of the 
dominant X-point (this distance is chosen to be outside 
of the initial field-reversal region of thickness A). Here, 



FIG. 14. Normalised reconnection rates E* against time for 
fluid Harris-sheet simulations with d e — 5 p s — 0.01, p = 10 -4 
and t]h = 10 -9 . Shown are L x /L x 0 = 0.5 (orange solid, same 
as in Fig. 10), L x /L x 0 = 2 (purple, dashed), and L x /L x 0 = 8 


(green, dotted) where L x 0 = 100d e 



Ex/E x o 


FIG. 15. Averaged peak rate < E* > against system-size 
L x /L x 0 , where L x 0 = lOO/i. From Harris-sheet fluid simula¬ 
tions with p — 10 -5 , t]h — 10 -9 , d e — 5 p s — 0.01 (orange A), 
p s = 5 d e = 0.01 (blue O)? an d fully kinetic simulations with 
d e = 5 p s = 0.01 (red ♦). Reprinted with permission from 
Stanier et ah, Phys. Plasmas, 22, 010701, (2015). Copyright 
(2015) American Institute of Physics. 


we average the rate (<>) over the time required for an 
Alfven wave to cross the whole box in the inflow direc¬ 
tion, so that the measurement of the peak rate is not 
overly affected by sharp spikes in the rate associated with 
secondary island formation, see Fig. [l4| Also, we neglect 
the final 25% of the simulation time when calculating this 
averaged peak rate as the upstream magnetic field can be¬ 
come significantly depleted close to saturation, leading to 
large increases in the normalised rate. 

In Fig [l5j we show results from: fluid simulations with 
p s = 0.01, d e = 0.002 (with FDWs); fluid runs with 
d e = 0.01, p s = 0.002 (without FDWs); and VPIC simu¬ 
lations with d e = 0.01 and p s = 0.002. The fluid simula¬ 
tions all have p = 10 -5 , t]h = 10 —9 , while the additional 
VPIC specific parameters are the same as in Sec. |V C[ 
but with l o P e/Qcex — 2.5 and v t h,e/c = 0.25. This change 
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in the VPIC specific parameters gives runs that are less 
expensive than the one in Sec. |V C[ particularly for the 
large simulations presented here, but the Alfven speed 
becomes mildly relativistic. Thus, we normalise by the 
relativistic Alfvenic rate for the VPIC runs, see Ref. |371 

Firstly, the case with FDWs has a clear system-size 
independent rate, with < E * 0.15. For the case 

without FDWs, there is some small variation in the peak 
rate that appears due to differences in secondary island 
formation between runs, even after the time averaging 
of the rate. However, there is no appreciable downward 
trend as the rate for the L x = 0.5 (total box length of 
100 d e ) simulation is very close to the L x = 8 (1600d e ) 
run. As the rate is a significant fraction of the Alfvenic 
rate for cases with and without FDWs, both cases sat¬ 
isfy all the criteria of fast reconnection; they are Alfvenic, 
independent of dissipation, and of system-size. In addi¬ 
tion, there is good agreement in the peak rates between 
VPIC and fluid calculations for the runs without FDWs. 
However, we note that for larger L x the time histories of 
do not agree quite as well as in Fig. 10, due to the 
sensitivity of secondary island formation on factors such 
as the level of noise (not shown). 

Previous low ~/3 kinetic simulations have found that the 
absolute values for the reconnection rates in simulations 
without FDWs can be as large asj^ or even exceed,^ 
those for simulations with FDWs. In this set of Harris- 
sheet simulations we find that the rate is typically larger 
for the simulations with p s = 5 d e = 0.01 (with FDWs) 
than those with d e = 5 p s = 0.01 (without FDWs), al¬ 
though this difference is within a factor of two. We also 
note that for the island coalescence runs in Fig. [4j and 
between the Harris sheet runs with A = 0.5 in Figs. p~8| i 
and [T8| ) of Sec. m the rates are much closer between 
the two cases. 


Figures [16] and [17) show several time snapshots of the 
current density and magnetic flux from PIXIE2D fluid 
simulations in cases without (top halves of each time 
snapshot) and with (bottom halves of each time snap¬ 
shot) FDWs, for the L x = 2 and L x = 8 runs respec¬ 
tively. For L x — 2 there are clear qualitative differences 
at t > 4, after the initial perturbation has compressed 
the current layer. For the case without FDWs the cur¬ 
rent layer shows multiple islands at t = 4, and a large 
island at t = 6,8, which appears to limit w v . In contrast 
to this, the current layer in the case with FDWs does not 
produce a secondary island until t « 8, where the recon¬ 
nection is close to saturating. The structure of the flux 
and current in the case with FDWs does resemble the 
open X-point configuration® more closely than the case 
without FDWs in this set of runs. Such an open X-point 
configuration has been linked td 44 * 62 ^ the conservation of 
the generalised fields G± — ^ — d^j ± d e p s uj , which in 
the collisionless limit (tjh = p = 0) are Lagrangian con¬ 
served quanities in the flow fields v± = z x V</4, with 
(j>± = 4> ± p s ^/d e . Here we have finite dissipation so 
that these generalised fields can reconnect, but we do see 
that G± is advected by the counter-twisted flows v± to 


give localised current and vorticity structures along the 
separatrices of <// ± , as in Ref. [44] Interestingly, Figure 
5 of Ref. 44] shows that the X-point begins to open up 
even when p s = d e and there are no FDWs present. The 
opening angle between the separatrices has been shown 
to increase in the presence of hot ions! 32! 

For the L x = 8 runs, the differences are less clear, as 
there are multiple islands at t = 8 in cases with and 
without FDWs. However, the opening angle between the 
separatrices at t = 8 does appear larger in the case with 
FDWs, whereas there is more secondary island formation 
at late times (t > 8) in the case without FDWs. Finally, 
despite these qualitative differences, the time at which 
reconnection saturates is very similar between cases with 
(t ~ 20) and without (t ~ 25) FDWs, at which point a 
similar amount of flux has been reconnected. 

It is not clear whether the mechanism for limiting w v 
to give system-size independent reconnection rates dif¬ 
fers between the cases with and without FDWs. There 
are significant qualitative differences in structure of the 
current layer between these cases in Figs. [l6j [I7| where 
it appears that secondary island formation plays a role 
in regulating the length of the layer in the case with¬ 
out FDWs. However, the reconnection rates in these 
L x = 2,8 runs without FDWs are within a factor of two 
of the rates for the runs with FDWs, and also similar 
to the L x = 0.5 fluid simulation without FDWs, where 
there is no secondary island formation. 


VII. COMPARISON WITH PREVIOUS TWO-FLUID 
STUDY 

In this section, the results from a set of two-fluid and 
kinetic simulations are presented that can be more closely 
compared with earlier two-fluid results from Fig. 3a of 
Ref. [27] Here, fluid simulations are performed in a 
quarter-box (x,y) G [0,10]x [0,5] using the Harris-sheet 
initial conditions of Eq. pi) ). Two different values for 
the equilibrium current sheet thickness are used, A = 0.5 
and A = 1. 

Figure 3a of Ref. 27 shows the reconnection rate 
against /3 X /2 (the plasma -(3 defined in terms of the up¬ 
stream reconnecting field B x 0 ) for a fixed guide to recon¬ 
necting field ratio of Bq/B x o = 30, and a fixed mass-ratio 
of rrii/m e = 82.6. Here, we match these parameters in 
our reduced model by choosing d e = 1 /\/rrii/m e = 0.11, 
and vary since p s oc (Z^/2) 1 / 2 with Bq/B x q and 
rrii/m e fixed. The threshold for the plasma to support 
fast-dispersive kinetic Alfven waves is p s > d e as before, 
which corresponds to /3 x /2 > (m e /rrii)(BQ /H 2 0 ) = 10.9, 
in agreement with the threshold (dotted line) given in 
Fig. 3a of Ref. El For the PIXIE2D fluid simulations, 
we use t]h = 10 -6 ,/i = 10 -4 , which ensure that the DR 
thickness is sufficiently below the two-fluid scales. 

Fig. [T8 } i shows the reconnection rate E z , measured at 
the dominant X-point and normalised by the Alfvenic 
rate at the upstream boundary H 2 0 , plotted against p s . 
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FIG. 16. Current and flux for L x — 2 from PIXIE2D fluid simulations. Top half of each snapshot is for d e — 0.01, p s = 0.002 
(no FDWs), and bottom half is p s = 0.01, d e = 0.002 (FDWs). 


Here, blue y are from PIXIE2D fluid simulations with 
initial current sheet thickness A = 0.5 (the thickness used 
in Ref. [27] h blue O are from PIXIE2D fluid simulations 
with A = 1, and green □ are from VPIC simulations with 
A = di = 1, Bq/B x q = 30, rrii/m e = 82.6 and varying 
p s . As was found in Ref. 27 ; we see the reconnection rate 
decrease with p s oc (/3 a ,/2) 1 / 2 such that the rate is lower 
in the region where there are no FDWs than with FDWs. 
However, there are some subtle differences between the 
A = 0.5 runs presented here compared with Ref. 127 that 
are worth highlighting. 

Firstly, it is mentioned in Ref. [27]that there is a factor 
of two decrease in the rate between a regime with both 
whistler and kinetic Alfven waves and a regime with ki¬ 
netic Alfven waves only. In Fig. 3a of Ref. l27l the max¬ 
imum P x /2 = 30 which corresponds to p s ~ 0.18 in our 
Fig. [18] a. We find that for the simulation with larger 


p s = 0.33, such that p s / A is closer to unity, the rate 
E z ^0.1 even without whistler waves. 

Secondly, we find that the decrease in the rate with p s 
in the region p s < d e (No FDWs) is slightly flatter than 
that found in Ref. im such that for p s = f3 x = 0 we find 
E z « 0.033, which is more than twice the value reported 
in Ref. 27. In these simulations (with A = 0.5) we do 
see some secondary islands for the p s = 0,0.022,0.044 
runs that may effect the rates in this region. Note that 
for these runs we repeat the simulation in the half-box 
(x,y) e [-10,10] x [0, 5] to prevent islands from being 
trapped at x = 0. 

We do not see any secondary island formation in any 
of the simulations reported here with the thicker initial 
current sheets A = 1 (blue O i n Fig. [l8]a). The peak 
rates for these runs are smaller in magnitude, but we 
see the same overall trend with a flattening of the rate 
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Max. 




Man. 53.94 
Min: -45^7 


FIG. 17. Current and flux for L x — 8 from PIXIE2D fluid simulations. Top half of each snapshot is for d e — 0.01, p s — 0.002 
(no FDWs), and bottom half is p s m 0.01, d e = 0.002 (FDWs). 


when p s < d e , which is also shown in fully kinetic VPIC 
simulations (green □), and there is a sharper increase in 
the rate when p s > d e . This trend is consistent with 
the idea that the reconnection rate depends only on the 
maximum of the two-fluid scales h = max[p s , d e \. 


Figure [l8]b shows the reconnection rate from a set of 
Harris sheet simulations with A = 1 using the same set¬ 
up, but with fixed p s = 0.11 (dashed line) and varying d e . 
For these simulations, kinetic Alfven waves are supported 
to the left of the dashed line, while there are no FDWs 
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Ps 



d e 

FIG. 18. a) Reconnection rate E z from Harris-sheet simula¬ 
tions with fixed d e = 0.11 and different p s . Shown are fluid 
simulations with initial current sheet thickness A = 0.5 (blue 
V), and fluid (blue O) an d kinetic (green □) simulations with 
A = 1. b) E z for fixed p s — 0.11 and varying d e for fluid (ma¬ 
genta filled □) and kinetic (green hollow □) simulations with 
A = 1. 


to the right of the dashed line. Remarkably, the peak 
rates increase with increasing d e , with almost the same 
values as for increasing p s with fixed d e in Fig. [18] a, 
supporting the conclusion that the rate depends only on 
the maximum of the two-fluid scales, regardless of the 
presence of FDWs. Fully kinetic (VPIC) runs with the 
same fixed p s and varying d e also show the same trend 
(green □). 


To understand why the rate increases for fixed p s and 
increasing d e , Fig. [T9] shows the two-fluid region param¬ 
eters that govern the rate in Eq. (15). Interestingly, 
Fig. 19 r shows the aspect-ratio of the outer region d e /w v , 
at the time of peak rate, remains roughly constant with 
value « 0.1 over a large range in d e . However, the 
Alfvenic rate on the edge of the two-fluid region at the 
time of peak rate, B^ hJ in Fig. 19 3 does increase with d e . 
This accounts for the scaling of the reconnection rate, 
as shown in 19", where the expression for the rate in 
Eq. (Il5|) has been scaled by a numerical factor of 3.5. 


We find that for runs with a smaller ratio of the max¬ 
imum two-fluid scale to equilibrium current sheet thick¬ 
ness, h/ A, there is a longer onset period before the peak 
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FIG. 19. a) Aspect ratio of outer two-fluid region d e /w v , b) 
Alfvenic rate, B^ hl on edge of two-fluid region, c) Actual rate 
at X-point E z (magenta □), and rate from Eq. ( |19| ) multiplied 
by 3.5 (blue O)- All results from A = 1 fluid simulations with 
fixed p s = 0.11. 


reconnection rate is reached, and the upstream Alfvenic 
rate B^ h has time to significantly decrease. This is very 
similar to the results found in a previous study of Hall- 
MHD reconnection with large equilibrium current sheet 
thickness.^ In both cases, with and without FDWs, if 
the peak rate is normalised by the instantaneous value of 
B^ h , then it remains roughly constant against h. From 
the results presented in this section, we argue that the 
presence or absence of FDWs does not play a primary 
role in whether the magnitude of the reconnection rate 
is large, as was concluded previously in Ref. [27l 
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VIII. SUMMARY AND CONCLUSIONS 


In this paper we show that low-/3 two-fluid reconnec¬ 
tion is formally fast, regardless of the presence of fast- 
dispersive waves (FDWs) that have previously been sug¬ 
gested to play a critical roleP^ This paper gives addi¬ 
tional details on the method used, and gives additional 
evidence to support the conclusions in a recently pub¬ 
lished letterWe show that reconnection in Harris sheet 
geometry is independent of the physics of the dissipation 
region even in the absence of FDWs, and demonstrate 
for cases with and without FDWs that the reconnection 
rate and the overall length of the current layer are in very 
good agreement between fluid and kinetic simulations. 

For simulations with very large system-sizes, it is not 
yet clear which mechanism limits the length of the layer 
w v and, in particular, whether the mechanism is the same 


regardless of the presence of FDWs. In Sec. VI we do see 
qualitative differences between the structure of the layer, 
such as a larger opening angle between the separatrices in 
the case with FDWs, and more secondary islands in the 
case without FDWs. However, the peak rates between 
the two cases agree to within a factor of two, and the 
time until saturation is also very similar. Furthermore, in 
the case without FDWs, the reconnection rates for small 
systems (L x = 0.5) with no secondary islands are the 
same as that for the larger systems (L x = 2,8) which do 
have secondary islands. Interestingly, for the Harris-sheet 
simulations without FDWs presented in Fig. [l9j there is 
no secondary island formation and the aspect ratio of the 
layer is found to be d e /w v « 0.1 across a large range in 
d e . This suggests that there is a mechanism other than 
secondary island formation that can limit the length of 
the layer w v in the case without FDWs. 

Also in Sec. m we argue that the decrease in the re¬ 
connection rate with f3 x /2 is not related to the presence 
of fast-dispersive waves, which was concluded in a previ¬ 
ously two-fluid study. 27 Instead, we argue that it is due 
to the decrease in the upstream Alfvenic rate B^ h asso¬ 
ciated with a longer onset period for runs with smaller 
values of h/X. This phenomenon is similar to that found 
in a previous Hall-MHD study, which considered how the 
reconnection rate scaled for small d{/ A. 63 
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